Identification of the potassium-binding site in serotonin transporter.

Clearance of serotonin (5-hydroxytryptamine, 5-HT) from the synaptic cleft after neuronal signaling is mediated by serotonin transporter (SERT), which couples this process to the movement of a Na+ ion down its chemical gradient. After release of 5-HT and Na+ into the cytoplasm, the transporter faces a rate-limiting challenge of resetting its conformation to be primed again for 5-HT and Na+ binding. Early studies of vesicles containing native SERT revealed that K+ gradients can provide an additional driving force, via K+ antiport. Moreover, under appropriate conditions, a H+ ion can replace K+. Intracellular K+ accelerates the resetting step. Structural studies of SERT have identified two binding sites for Na+ ions, but the K+ site remains enigmatic. Here, we show that K+ antiport can drive substrate accumulation into vesicles containing SERT extracted from a heterologous expression system, allowing us to study the residues responsible for K+ binding. To identify candidate binding residues, we examine many cation binding configurations using molecular dynamics simulations, predicting that K+ binds to the so-called Na2 site. Site-directed mutagenesis of residues in this site can eliminate the ability of both K+ and H+ to drive 5-HT accumulation into vesicles and, in patch clamp recordings, prevent the acceleration of turnover rates and the formation of a channel-like state by K+ or H+. In conclusion, the Na2 site plays a pivotal role in orchestrating the sequential binding of Na+ and then K+ (or H+) ions to facilitate 5-HT uptake in SERT.

VMAT structures reveal exciting targets for drug development.

Vesicular monoamine transporter (VMAT)-2 has a crucial role in the neurotransmission of biogenic amines. Recently, Dalton et al., Pidathala et al., Wu et al., and Wang et al. individually reported cryo-electron microscopy (EM) structures of human VMAT2, offering opportunities for developing improved therapeutics and deep insights into the functioning of this protein.

EncoMPASS: An encyclopedia of membrane proteins analyzed by structure and symmetry.

Protein structure determination and prediction, active site detection, and protein sequence alignment techniques all exploit information about protein structure and structural relationships. For membrane proteins, however, there is limited agreement among available online tools for highlighting and mapping such structural similarities. Moreover, no available resource provides a systematic overview of quaternary and internal symmetries, and their orientation relative to the membrane, despite the fact that these properties can provide key insights into membrane protein function and evolution. Here, we describe the Encyclopedia of Membrane Proteins Analyzed by Structure and Symmetry (EncoMPASS), a database for relating integral membrane proteins of known structure from the points of view of sequence, structure, and symmetry. EncoMPASS is accessible through a web interface, and its contents can be easily downloaded. This allows the user not only to focus on specific proteins, but also to study general properties of the structure and evolution of membrane proteins.

Membrane free-energy landscapes derived from atomistic dynamics explain nonuniversal cholesterol-induced stiffening.

All lipid membranes have inherent morphological preferences and resist deformation. Yet adaptations in membrane shape can and do occur at multiple length scales. While this plasticity is crucial for cellular physiology, the factors controlling the morphological energetics of lipid bilayers and the dominant mechanisms of membrane remodeling remain to be fully understood. An ongoing debate regarding the universality of the stiffening effect of cholesterol underscores the challenges facing this field, both experimentally and theoretically, even for simple lipid mixtures. On the computational side, we have argued that enhanced-sampling all-atom molecular dynamics simulations are uniquely suited for the quantification of membrane conformational energetics, as they minimize a priori assumptions and permit analysis of bilayers in deformed states. To showcase this approach, we examine reported inconsistencies between alternative experimental measurements of bending moduli for cholesterol-enriched membranes. Specifically, we analyze lipid bilayers with different chain saturation and compute free-energy landscapes for curvature deformations distributed over areas from ∼5 to ∼60 nm2. These enhanced simulations, totaling over 100 µs of sampling time, enable us to directly quantify both bending and tilt moduli and to dissect the contributing factors and molecular mechanisms of curvature generation at each length scale. Our results show that the effects of cholesterol on bending rigidity are lipid-specific and suggest that this specificity arises from differences in the torsional dynamics of the acyl chains. In summary, we demonstrate that quantitative relationships can now be established between lipid structure and bending energetics, paving the way for addressing open fundamental questions in cell membrane mechanics.

Membrane free-energy landscapes derived from atomistic dynamics explain nonuniversal cholesterol-induced stiffening.

All lipid membranes have inherent morphological preferences and resist deformation. Yet adaptations in membrane shape can and do occur at multiple length scales. While this plasticity is crucial for cellular physiology, the factors controlling the morphological energetics of lipid bilayers and the dominant mechanisms of membrane remodeling remain unclear. An ongoing debate regarding the universality of the stiffening effect of cholesterol underscores the challenges facing this field, both experimentally and theoretically, even for simple lipid mixtures. On the computational side, we have argued that enhanced- sampling all-atom molecular dynamics simulations are uniquely suited for quantification of membrane conformational energetics, not only because they minimize a-priori assumptions, but also because they permit analysis of bilayers in deformed states. To showcase this approach, we examine reported inconsistencies between alternative experimental measurements of bending moduli for cholesterol-enriched membranes. Specifically, we analyze lipid bilayers with different chain saturation, and compute free-energy landscapes for curvature deformations distributed over areas from ∼5 to ∼60 nm 2 . These enhanced simulations, totaling over 100 microseconds of sampling time, enable us to directly quantify both bending and tilt moduli, and to dissect the contributing factors and molecular mechanisms of curvature generation at each length scale. Our results show that cholesterol effects are lipid-specific, in agreement with giantvesicle measurements, and explain why experiments probing nanometer scale lipid dynamics diverge. In summary, we demonstrate that quantitative structure-mechanics relationships can now be established for heterogenous membranes, paving the way for addressing open fundamental questions in cell membrane mechanics. Significance: Elucidating the energetics and mechanisms of membrane remodeling is an essential step towards understanding cell physiology. This problem is challenging, however, because membrane bending involves both large-scale and atomic-level dynamics, which are difficult to measure simultaneously. A recent controversy regarding the stiffening effect of cholesterol, which is ubiquitous in animal cells, illustrates this challenge. We show how enhanced molecular-dynamics simulations can bridge this length-scale gap and reconcile seemingly incongruent observations. This approach facilitates a conceptual connection between lipid chemistry and membrane mechanics, thereby providing a solid basis for future research on remodeling phenomena, such as in membrane trafficking or viral infection.

Insights into autoregulation of a membrane protein complex by its cytoplasmic domains.

Membrane transporters mediate the passage of molecules across membranes and are essential for cellular function. While the transmembrane region of these proteins is responsible for substrate transport, often the cytoplasmic regions are required for modulating their activity. However, it can be difficult to obtain atomic-resolution descriptions of these autoregulatory domains by classical structural biology techniques, especially if they lack a single, defined structure. The betaine permease, BetP, a homotrimer, is a prominent and well-studied example of a membrane protein whose autoregulation depends on cytoplasmic N- and C-terminal segments. These domains sense and transduce changes in K+ concentration and in lipid bilayer properties caused by osmotic stress. However, structural data for these terminal domains is incomplete, which hinders a clear description of the molecular mechanism of autoregulation. Here we used microsecond-scale molecular simulations of the BetP trimer to compare reported conformations of the 45-amino-acid long C-terminal tails. The simulations provide support for the idea that the conformation derived from electron microscopy (EM) data represents a more stable global orientation of the C-terminal segment under downregulating conditions while also providing a detailed molecular description of its dynamics and highlighting specific interactions with lipids, ions, and neighboring transporter subunits. A missing piece of the molecular puzzle is the N-terminal segment, whose dynamic nature has prevented structural characterization. Using Rosetta to generate ensembles of de novo conformations in the context of the EM-derived structure robustly identifies two features of the N-terminal tail, namely 1) short helical elements and 2) an orientation that would confine potential interactions to the protomer in the counterclockwise direction (viewed from the cytoplasm). Since each C-terminal tail only contacts the protomer in the clockwise direction, these results indicate an intricate interplay between the three protomers of BetP in the downregulated protein and a multidirectionality that may facilitate autoregulation of transport.

AlignMe: an update of the web server for alignment of membrane protein sequences.

The AlignMe web server is dedicated to accurately aligning sequences of membrane proteins, a particularly challenging task due to the strong evolutionary divergence and the low compositional complexity of hydrophobic membrane-spanning proteins. AlignMe can create pairwise alignments of either two primary amino acid sequences or two hydropathy profiles. The web server for AlignMe has been continuously available for >10 years, supporting 1000s of users per year. Recent improvements include anchoring, multiple submissions, and structure visualization. Anchoring is the ability to constrain a position in an alignment, which allows expert information about related residues in proteins to be incorporated into an alignment without manual modification. The original web interface to the server limited the user to one alignment per submission, hindering larger scale studies. Now, batches of alignments can be initiated with a single submission. Finally, to provide structural context for the relationship between proteins, sequence similarity can now be mapped onto one or more structures (or structural models) of the proteins being aligned, by links to MutationExplorer, a web-based visualization tool. Together with a refreshed user interface, these features further enhance an important resource in the membrane protein community. The AlignMe web server is freely available at https://www.bioinfo.mpg.de/AlignMe/.

On the Role of a Conserved Methionine in the Na+-Coupling Mechanism of a Neurotransmitter Transporter Homolog.

Excitatory amino acid transporters (EAAT) play a key role in glutamatergic synaptic communication. Driven by transmembrane cation gradients, these transporters catalyze the reuptake of glutamate from the synaptic cleft once this neurotransmitter has been utilized for signaling. Two decades ago, pioneering studies in the Kanner lab identified a conserved methionine within the transmembrane domain as key for substrate turnover rate and specificity; later structural work, particularly for the prokaryotic homologs GltPh and GltTk, revealed that this methionine is involved in the coordination of one of the three Na+ ions that are co-transported with the substrate. Albeit extremely atypical, the existence of this interaction is consistent with biophysical analyses of GltPh showing that mutations of this methionine diminish the binding cooperativity between substrates and Na+. It has been unclear, however, whether this intriguing methionine influences the thermodynamics of the transport reaction, i.e., its substrate:ion stoichiometry, or whether it simply fosters a specific kinetics in the binding reaction, which, while influential for the turnover rate, do not fundamentally explain the ion-coupling mechanism of this class of transporters. Here, studies of GltTk using experimental and computational methods independently arrive at the conclusion that the latter hypothesis is the most plausible, and lay the groundwork for future efforts to uncover the underlying mechanism.

Modeling the native ensemble of PhuS using enhanced sampling MD and HDX-ensemble reweighting.

The cytoplasmic heme binding protein from Pseudomonas aeruginosa, PhuS, plays two essential roles in regulating heme uptake and iron homeostasis. First, PhuS shuttles exogenous heme to heme oxygenase (HemO) for degradation and iron release. Second, PhuS binds DNA and modulates the transcription of the prrF/H small RNAs (sRNAs) involved in the iron-sparing response. Heme binding to PhuS regulates this dual function, as the unliganded form binds DNA, whereas the heme-bound form binds HemO. Crystallographic studies revealed nearly identical structures for apo- and holo-PhuS, and yet numerous solution-based measurements indicate that heme binding is accompanied by large conformational rearrangements. In particular, hydrogen-deuterium exchange mass spectrometry (HDX-MS) of apo- versus holo-PhuS revealed large differences in deuterium uptake, notably in α-helices 6, 7, and 8 (α6,7,8), which contribute to the heme binding pocket. These helices were mostly labile in apo-PhuS but largely protected in holo-PhuS. In contrast, in silico-predicted deuterium uptake levels of α6,7,8 from molecular dynamics (MD) simulations of the apo- and holo-PhuS structures are highly similar, consistent only with the holo-PhuS HDX-MS data. To rationalize this discrepancy between crystal structures, simulations, and observed HDX-MS, we exploit a recently developed computational approach (HDXer) that fits the relative weights of conformational populations within an ensemble of structures to conform to a target set of HDX-MS data. Here, a combination of enhanced sampling MD, HDXer, and dimensionality reduction analysis reveals an apo-PhuS conformational landscape in which α6, 7, and 8 are significantly rearranged compared to the crystal structure, including a loss of secondary structure in α6 and the displacement of α7 toward the HemO binding interface. Circular dichroism analysis confirms the loss of secondary structure, and the extracted ensembles of apo-PhuS and of heme-transfer-impaired H212R mutant, are consistent with known heme binding and transfer properties. The proposed conformational landscape provides structural insights into the modulation by heme of the dual function of PhuS.

Refining pairwise sequence alignments of membrane proteins by the incorporation of anchors.

The alignment of primary sequences is a fundamental step in the analysis of protein structure, function, and evolution, and in the generation of homology-based models. Integral membrane proteins pose a significant challenge for such sequence alignment approaches, because their evolutionary relationships can be very remote, and because a high content of hydrophobic amino acids reduces their complexity. Frequently, biochemical or biophysical data is available that informs the optimum alignment, for example, indicating specific positions that share common functional or structural roles. Currently, if those positions are not correctly matched by a standard pairwise sequence alignment procedure, the incorporation of such information into the alignment is typically addressed in an ad hoc manner, with manual adjustments. However, such modifications are problematic because they reduce the robustness and reproducibility of the aligned regions either side of the newly matched positions. Previous studies have introduced restraints as a means to impose the matching of positions during sequence alignments, originally in the context of genome assembly. Here we introduce position restraints, or "anchors" as a feature in our alignment tool AlignMe, providing an aid to pairwise global sequence alignment of alpha-helical membrane proteins. Applying this approach to realistic scenarios involving distantly-related and low complexity sequences, we illustrate how the addition of anchors can be used to modify alignments, while still maintaining the reproducibility and rigor of the rest of the alignment. Anchored alignments can be generated using the online version of AlignMe available at www.bioinfo.mpg.de/AlignMe/.

Chloride-dependent conformational changes in the GlyT1 glycine transporter.

The human GlyT1 glycine transporter requires chloride for its function. However, the mechanism by which Cl- exerts its influence is unknown. To examine the role that Cl- plays in the transport cycle, we measured the effect of Cl- on both glycine binding and conformational changes. The ability of glycine to displace the high-affinity radioligand [3H]CHIBA-3007 required Na+ and was potentiated over 1,000-fold by Cl- We generated GlyT1b mutants containing reactive cysteine residues in either the extracellular or cytoplasmic permeation pathways and measured changes in the reactivity of those cysteine residues as indicators of conformational changes in response to ions and substrate. Na+ increased accessibility in the extracellular pathway and decreased it in the cytoplasmic pathway, consistent with stabilizing an outward-open conformation as observed in other members of this transporter family. In the presence of Na+, both glycine and Cl- independently shifted the conformation of GlyT1b toward an outward-closed conformation. Together, Na+, glycine, and Cl- stabilized an inward-open conformation of GlyT1b. We then examined whether Cl- acts by interacting with a conserved glutamine to allow formation of an ion pair that stabilizes the closed state of the extracellular pathway. Molecular dynamics simulations of a GlyT1 homolog indicated that this ion pair is formed more frequently as that pathway closes. Mutation of the glutamine blocked the effect of Cl-, and substituting it with glutamate or lysine resulted in outward- or inward-facing transporter conformations, respectively. These results provide an unexpected insight into the role of Cl- in this family of transporters.

Interpreting Hydrogen-Deuterium Exchange Experiments with Molecular Simulations: Tutorials and Applications of the HDXer Ensemble Reweighting Software [Article v1.0].

Hydrogen-deuterium exchange (HDX) is a comprehensive yet detailed probe of protein structure and dynamics and, coupled to mass spectrometry, has become a powerful tool for investigating an increasingly large array of systems. Computer simulations are often used to help rationalize experimental observations of exchange, but interpretations have frequently been limited to simple, subjective correlations between microscopic dynamical fluctuations and the observed macroscopic exchange behavior. With this in mind, we previously developed the HDX ensemble reweighting approach and associated software, HDXer, to aid the objective interpretation of HDX data using molecular simulations. HDXer has two main functions; first, to compute H-D exchange rates that describe each structure in a candidate ensemble of protein structures, for example from molecular simulations, and second, to objectively reweight the conformational populations present in a candidate ensemble to conform to experimental exchange data. In this article, we first describe the HDXer approach, theory, and implementation. We then guide users through a suite of tutorials that demonstrate the practical aspects of preparing experimental data, computing HDX levels from molecular simulations, and performing ensemble reweighting analyses. Finally we provide a practical discussion of the capabilities and limitations of the HDXer methods including recommendations for a user's own analyses. Overall, this article is intended to provide an up-to-date, pedagogical counterpart to the software, which is freely available at https://github.com/Lucy-Forrest-Lab/HDXer.

Global alignment and assessment of TRP channel transmembrane domain structures to explore functional mechanisms.

The recent proliferation of published TRP channel structures provides a foundation for understanding the diverse functional properties of this important family of ion channel proteins. To facilitate mechanistic investigations, we constructed a structure-based alignment of the transmembrane domains of 120 TRP channel structures. Comparison of structures determined in the absence or presence of activating stimuli reveals similar constrictions in the central ion permeation pathway near the intracellular end of the S6 helices, pointing to a conserved cytoplasmic gate and suggesting that most available structures represent non-conducting states. Comparison of the ion selectivity filters toward the extracellular end of the pore supports existing hypotheses for mechanisms of ion selectivity. Also conserved to varying extents are hot spots for interactions with hydrophobic ligands, lipids and ions, as well as discrete alterations in helix conformations. This analysis therefore provides a framework for investigating the structural basis of TRP channel gating mechanisms and pharmacology, and, despite the large number of structures included, reveals the need for additional structural data and for more functional studies to establish the mechanistic basis of TRP channel function.

Interpretation of HDX Data by Maximum-Entropy Reweighting of Simulated Structural Ensembles.

Hydrogen-deuterium exchange combined with mass spectrometry (HDX-MS) is a widely applied biophysical technique that probes the structure and dynamics of biomolecules without the need for site-directed modifications or bio-orthogonal labels. The mechanistic interpretation of HDX data, however, is often qualitative and subjective, owing to a lack of quantitative methods to rigorously translate observed deuteration levels into atomistic structural information. To help address this problem, we have developed a methodology to generate structural ensembles that faithfully reproduce HDX-MS measurements. In this approach, an ensemble of protein conformations is first generated, typically using molecular dynamics simulations. A maximum-entropy bias is then applied post hoc to the resulting ensemble such that averaged peptide-deuteration levels, as predicted by an empirical model, agree with target values within a given level of uncertainty. We evaluate this approach, referred to as HDX ensemble reweighting (HDXer), for artificial target data reflecting the two major conformational states of a binding protein. We demonstrate that the information provided by HDX-MS experiments and by the model of exchange are sufficient to recover correctly weighted structural ensembles from simulations, even when the relevant conformations are rarely observed. Degrading the information content of the target data-e.g., by reducing sequence coverage, by averaging exchange levels over longer peptide segments, or by incorporating different sources of uncertainty-reduces the structural accuracy of the reweighted ensemble but still allows for useful insights into the distinctive structural features reflected by the target data. Finally, we describe a quantitative metric to rank candidate structural ensembles according to their correspondence with target data and illustrate the use of HDXer to describe changes in the conformational ensemble of the membrane protein LeuT. In summary, HDXer is designed to facilitate objective structural interpretations of HDX-MS data and to inform experimental approaches and further developments of theoretical exchange models.

MemSTATS: A Benchmark Set of Membrane Protein Symmetries and Pseudosymmetries.

In membrane proteins, symmetry and pseudosymmetry often have functional or evolutionary implications. However, available symmetry detection methods have not been tested systematically on this class of proteins because of the lack of an appropriate benchmark set. Here we present MemSTATS, a publicly available benchmark set of both quaternary- and internal-symmetries in membrane protein structures. The symmetries are described in terms of order, repeated elements, and orientation of the axis with respect to the membrane plane. Moreover, using MemSTATS, we compare the performance of four widely used symmetry detection algorithms and highlight specific challenges and areas for improvement in the future.

Large-scale state-dependent membrane remodeling by a transporter protein.

That channels and transporters can influence the membrane morphology is increasingly recognized. Less appreciated is that the extent and free-energy cost of these deformations likely varies among different functional states of a protein, and thus, that they might contribute significantly to defining its mechanism. We consider the trimeric Na+-aspartate symporter GltPh, a homolog of an important class of neurotransmitter transporters, whose mechanism entails one of the most drastic structural changes known. Molecular simulations indicate that when the protomers become inward-facing, they cause deep, long-ranged, and yet mutually-independent membrane deformations. Using a novel simulation methodology, we estimate that the free-energy cost of this membrane perturbation is in the order of 6-7 kcal/mol per protomer. Compensating free-energy contributions within the protein or its environment must thus stabilize this inward-facing conformation for the transporter to function. We discuss these striking results in the context of existing experimental observations for this and other transporters.

RNA Granules Hitchhike on Lysosomes for Long-Distance Transport, Using Annexin A11 as a Molecular Tether.

Long-distance RNA transport enables local protein synthesis at metabolically-active sites distant from the nucleus. This process ensures an appropriate spatial organization of proteins, vital to polarized cells such as neurons. Here, we present a mechanism for RNA transport in which RNA granules "hitchhike" on moving lysosomes. In vitro biophysical modeling, live-cell microscopy, and unbiased proximity labeling proteomics reveal that annexin A11 (ANXA11), an RNA granule-associated phosphoinositide-binding protein, acts as a molecular tether between RNA granules and lysosomes. ANXA11 possesses an N-terminal low complexity domain, facilitating its phase separation into membraneless RNA granules, and a C-terminal membrane binding domain, enabling interactions with lysosomes. RNA granule transport requires ANXA11, and amyotrophic lateral sclerosis (ALS)-associated mutations in ANXA11 impair RNA granule transport by disrupting their interactions with lysosomes. Thus, ANXA11 mediates neuronal RNA transport by tethering RNA granules to actively-transported lysosomes, performing a critical cellular function that is disrupted in ALS.

Intestinal serotonin and fluoxetine exposure modulate bacterial colonization in the gut.

The gut microbiota regulates levels of serotonin (5-hydroxytryptamine (5-HT)) in the intestinal epithelium and lumen1-5. However, whether 5-HT plays a functional role in bacteria from the gut microbiota remains unknown. We demonstrate that elevating levels of intestinal lumenal 5-HT by oral supplementation or genetic deficiency in the host 5-HT transporter (SERT) increases the relative abundance of spore-forming members of the gut microbiota, which were previously reported to promote host 5-HT biosynthesis. Within this microbial community, we identify Turicibacter sanguinis as a gut bacterium that expresses a neurotransmitter sodium symporter-related protein with sequence and structural homology to mammalian SERT. T. sanguinis imports 5-HT through a mechanism that is inhibited by the selective 5-HT reuptake inhibitor fluoxetine. 5-HT reduces the expression of sporulation factors and membrane transporters in T. sanguinis, which is reversed by fluoxetine exposure. Treating T. sanguinis with 5-HT or fluoxetine modulates its competitive colonization in the gastrointestinal tract of antibiotic-treated mice. In addition, fluoxetine reduces the membership of T. sanguinis in the gut microbiota of conventionally colonized mice. Host association with T. sanguinis alters intestinal expression of multiple gene pathways, including those important for lipid and steroid metabolism, with corresponding reductions in host systemic triglyceride levels and inguinal adipocyte size. Together, these findings support the notion that select bacteria indigenous to the gut microbiota signal bidirectionally with the host serotonergic system to promote their fitness in the intestine.

A structural model of the human serotonin transporter in an outward-occluded state.

The human serotonin transporter hSERT facilitates the reuptake of its endogenous substrate serotonin from the synaptic cleft into presynaptic neurons after signaling. Reuptake regulates the availability of this neurotransmitter and therefore hSERT plays an important role in balancing human mood conditions. In 2016, the first 3D structures of this membrane transporter were reported in an inhibitor-bound, outward-open conformation. These structures revealed valuable information about interactions of hSERT with antidepressant drugs. Nevertheless, the question remains how serotonin facilitates the specific conformational changes that open and close pathways from the synapse and to the cytoplasm as required for transport. Here, we present a serotonin-bound homology model of hSERT in an outward-occluded state, a key intermediate in the physiological cycle, in which the interactions with the substrate are likely to be optimal. Our approach uses two template structures and includes careful refinement and comprehensive computational validation. According to microsecond-long molecular dynamics simulations, this model exhibits interactions between the gating residues in the extracellular pathway, and these interactions differ from those in an outward-open conformation of hSERT bound to serotonin. Moreover, we predict several features of this state by monitoring the intracellular gating residues, the extent of hydration, and, most importantly, protein-ligand interactions in the central binding site. The results illustrate common and distinct characteristics of these two transporter states and provide a starting point for future investigations of the transport mechanism in hSERT.